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Abstract. The steady current-voltage characteristics of ferromagnet-silicon- 
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in the presence of spin blockade generalizing the model developed by Pershin Yu V 
and Di Ventra M (2008 Phys. Rev. B 77 073301). The dependence of the critical 
current on both collector and emitter polarizations is obtained analytically. It is found 
that the current amplitude in the channel can be effectively tuned by varying the dif- 
ference between the collector and emitter ferromagnet polarizations which allows to 
perform the magnetic manipulation of the electrical current in wide class of both n- 
and p-doped, low- and high-Ohmic semiconductor channels coupled to ferromagnetic 
leads. 
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1. Introduction 

The progress of spintronics and physics of heterostructures which can be observed 
during the last years [U |2] is focused on various physical phenomena, and one of 
them which attracts a considerable attention is the spin-dependent transport through 
semiconductor/spin-polarized junctions PIlElEllTllHlElIiniiniinillSllIllIlS]. The 
physics of carrier polarization and its influence on transport in composite structures 
such as semi conductor /ferromagnet has been studied both theoretically [H El d El El 
[T2I [T3l [m |15] and experimentally [3l El |10l [11]. One of the models describing the 
spin-resolved carrier concentrations and currents at the junction is the two-component 
drift-diffusion model [6l [131 113] which predicted highly nonlinear and saturating current- 
voltage dependence at a single semi conductor /feromagnet junction due to the effect of 
spin blockade [131 E!- this model the detailed structure of the charge and current 
distribution at the junction area [HE] as well as the Schottky barriers [3 [8], the charge 
redistribution effects [9], and the bound states [T2] are not taken into consideration. 
Still, the quahtative and distinguishable behaviour of current saturation due to the 
effect of spin blockade is reliably predicted under various system parameters such as 
the junction/semiconductor resistance ratio. The spin blockade regime arises from the 
spatial distribution of the spin-minority carriers which cannot enter the ferromagnet 
region and form a cloud near the junction which growth prevents the further increase of 
spin-majority carrier transport if the current exceeds a threshold value called the critical 
current. Further studies have shown the importance and promising applications of this 
effect also for non- stationary phenomena such as spin memory effects [15]. The models 
described above were applied mainly to GaAs-based semiconductor channels, but is is 
known that the silicon-based structures are also of big interest for spintronics due to 
the dominating place of silicon in currently available electronic technologies. More, the 
technologies of fabricating the silicon/ferro magnet structures such as Si/Si:Mn formed 
on a basis of diluted magnetic semiconductors have been intensively developed during 
the last few years [Ml |T7j which makes their future applications in spintronics promising 
and creates certain questions about the phenomena described above. Is there a spin 
blockade regime in a silicon/ferromagnet junction at specific values of applied voltage, 
carrier mobility and concentration ? If so, what is the critical current density and how 
it depends on the silicon and ferromagnet parameters such as the carrier polarization 
in ferromagnets and the conductivity of the semiconductor channel? How deep can 
we modulate the current in the channel by manipulating the polarization of emitter 
or collector ferromagnets relative to each other? In the present manuscript we study 
these problems in the framework of a simple but effective model of transport in the spin 
blockade regime [131 [T3] which we generalize for the case of arbitrary carrier polarizations 
in the emitting and collecting ferromagnetic regions of the channel as well as for wide 
range of low- and high-Ohmic n-doped and p-doped silicon samples. It is found that the 
current can be deeply modulated by changing the spin alignment in the emitter and/or 
collector ferromagnet since the critical current density is very sensitive to it. We find the 



Current-voltage characteristics of tunable F-Si-F channels 



3 



analytical expression for the critical current density and calculate the current-voltage 
dependencies for various combinations of the channel/ contact resistance ratios, as well 
as for n- and p-type of doping with both high and low concentrations. The manuscript 
is organized as follows: in Section 2 we derive a model generalizing the description of 
the spin blockade regime for the two-ferromagnet channel with arbitrary polarizations in 
the emitter and collector ferromagnets and discuss the properties of the critical current 
density, in Section 3 we plot and discuss the current-voltage characteristics for various 
combinations of system parameters, and the conclusions are given in Section 4. 

2. The model 

The schematic view of the ferromagnet-silicon-ferromagnet channel is shown in Figure [TJ 
The collector ferromagnet with the junction resistance vi is separated from the emitter 
ferromagnet by a bulk silicon channel with length L which we consider as exceeding the 
spin diffusion length Ig given by [HI [131 E] 



where D and are the carrier diffusion coefficient and drift mobility, respectively, E 
is the electric field inside the channel, and r ~ 10 ns is the typical spin relaxation time 
[2]. One can see from ([T]) that Is is maximal at zero electric field when Ig = y/Drg 
and when the diffusion coefficient and spin relaxation time are big, or, since D is 
coupled to the mobility via the Einstein relation D = fiksT/e [18], it is clear that 
the spin diffusion length grows with the mobility. Since the mobility in the bulk silicon 
is typically lower than the one for GaAs [19] , for the given parameters of the silicon the 
condition L >> Igis already fulfilled if L > 10 mkm which is a reasonable channel length 
of bulk semiconductor structures. Hence, from the point view where the spin-resolved 
concentration decay length is considered, the distance between emitter and collector 
ferromagnets here can be taken as infinite which simplifies the boundary conditions. 
The key parameters of our model are the variable polarization degrees a and f5 in the 
collector and emitter ferromagnets describing the state of their non-ideality as well as 
the chosen direction of polarization for the majority of carriers. For example, the pair 
a = 1, P = 1/2 corresponds to the previously investigated case |14] where the carriers 
are fully polarized in the collector and fully unpolarized in the emitter. Correspondingly, 
the choice a = 0.8, (3 = 0.7 describes the situation when the collector and emitter are 
highly polarized in the same direction while the values a = 0.8, (3 = 0.3 describe the 
opposite polarization at the emitter. We adopt the labeling ji(2) for the spin-resolved 
current with spins in ferromagnets aligned up (down). If j is the total current density, 
than the spin-resolved current at the collector ji = aj which means that a = 1(0) 
corresponds to an ideal ferromagnet where magnetic moments are all aligned up (down), 
and a = 1/2 describes the unpolarized current. The same labeling is adopted for the 
carrier polarization generated by the emitter ferromagnet with polarization f3. The main 
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Figure 1. Schematic view of a ferromagnetic (F) - silicon (Si) - ferromagnetic 
(F) channel with variable polarization degrees a and (3 in collector and emitter 
ferromagnets. The channel and contact zero-bias resistances are tq and ri, and the 
direction of the carrier motion / is shown (for the electrons the actual current is — /). 
The spin blockade may occur in the contact region at the collector (left F) if the 
polarizations a and (3 significantly differ from each other, which enhances the contact 
resistance from the zero-bias value ri as long as the current builds up. 



purpose of the present paper is to find out how deep is the infiuence of the various a 
and (3 differences reflecting the polarization switch on the total electrical current in the 
channel. We include in our model the channel and contact zero-bias resistances r^ and vi 
(see Figure [1]) with variable ratio, and the direction of the carrier motion is labeled as / 
(for the electrons the actual current is — /). The spin blockade may occur in the contact 
region at the collector (left F in Figure [1]) if the polarizations a and (3 signiflcantly differ 
from each other, which enhances the contact resistance from the zero-bias value ri as 
long as the current builds up. 

The model for the spin-resolved current densities ji^2 and spin-resolved 
concentrations Ui + n2 = N where N is the total carrier concentration is well-known 
(HI [131 E! and consists of the continuity equation 

= ciivji,2 + ^("-2,1 - ^1,2) (2) 
and the equation for the spin-resolved current 

ji,2 = eni,2/U-E' + eDVni^2- (3) 

These equations should be accompanied by the boundary conditions at the collector 
contact X = and at the other boundary of the channel a^max = L » Ig which is in our 
problem can be considered as x = 00. In this case the steady-state solution of ([2]) exists 
[T3| [H] which has the form of two exponents decaying to their concentrations deflned 
by the boundary conditions at the emitter. In our model the emitter ferromagnet is 
described by the arbitrary polarization (3 which generalizes the unpolarized case /? = 1/2 
considered previously [131 IE] > so 



ni(oo) = (3N (4) 
712(00) = (l-/5)iV. (5) 

The steady-state solution of ([2]) satisfying the boundary condition (jS}) at the emitter 
and the normalizing condition Ui + n2 = N has the form 
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ni(x) =/3N-Ae-^'' 
n2(x) = (1 - /?)iV + Ae-^", 



(6) 
(7) 



where A = 1//^ is the inverse decay length defined in ([T]). The parameter A should 
be determined from the boundary condition at the collector when x = 0. If the collector 
is a non-ideal ferromagnet with the polarization a, then the spin-resolved current at the 
collector 



ji(0) = aj, 



(8) 



and the other boundary condition at the emitter ji(oo) = Pj is satisfied 
automatically. The total current density j is related to the electric field in the channel 
via usual relation j = efiNE which allows to construct a closed equation for the current 
density and applied voltage. From ([8]) it follows that the conductivity of the collector 
junction is proportional to the concentration of the majority carriers, and at the absence 
of the current the junction resistance has a predetermined value vi while the silicon 
channel with length L and cross-section S is described by the resistance defined in a usual 
way by the carrier mobility, concentration and geometric dimensions as tq = L/ (enNS). 
After a simple algebra we obtain the Ohmic law /(ro + ri(/)) = V which reads in our 
case similar to the one obtained previously [1A\ for a non-polarized emitter with f3 = 1/2, 



ro + ri 



1 



1 



V 



(9) 



but in our generalized model the parameter jc defined as the critical current density 
depends on both the collector and emitter ferromagnet polarizations a and jS: 



Jc Jc 



\( 



Q-/3/2 V 
/3 J 



(10) 



where 




(11) 



is the critical current density for the case of fully polarized collector ferromagnet 
with a = 1 and fully unpolarized emitter with (3 = 1/2 [T3]. Indeed, for a = 1 and 
(3 = 1/2 the value of jc from (fTOl) is equal to jjP^ . 

The main property of equation determining the current-voltage characteristic 
is the current saturation occurring at current densities comparable to jc- Thus, it is of 
interest to investigate the critical current density flTUl) in more details with respect to the 
tunable polarization parameters a and (3 of the collector and emitter ferromagnets. In 
Figure[2]the three-dimensional plot jc = jc{<^, P) is shown for the critical current density 
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Figure 2. (Colour online) Dependence of critical current density jc measured in units 
of jc"^ on the polarizations a and P in collector and emitter ferromagnets. The case of 
fully polarized carriers in the collector with a — 1 and fully unpolarized electrons in the 
emitter (3 = 1/2 considered in |14j is shown by a reference point A when jc = jc^"^ ■ The 
magnitude of jc diverges along the line a = f3 which means that in the case of a perfect 
coincidence of the spin alignment in the emitter and collector the critical current can 
never be reached and there is no spin blockade in this limit. In the opposite limit of 
maximum difference between a and /3 the value of jc is considerably smaller than ji*^^ , 
so the spin blockade here occurs at lower current densities. 

f fTOj) measured in units of jjP^ . The case of fully polarized carriers in the collector with 
a = 1 and fully unpolarized electrons in the emitter (3 = 1/2 considered in [T3] is shown 
by a reference point A when jc = jjP^. It is clear from (fTOj) that jc diverges along the 
line a = (3 which means that in the case of a perfect coincidence of the spin alignment 
in the emitter and collector the critical current can never be reached and there is no 
spin blockade in this limit. In the opposite limit of maximum difference between a and 
P the value of jc is considerably smaller than j^^\ so the spin blockade here occurs at 
lower current densities. The observed large variations of the critical current density are 
important for the desired high degree of control on the electrical current flowing through 
the ferromagnet-silicon-ferromagnet channel which can be achieved by variations of at 
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Figure 3. Current-voltage characteristics for a F-Si-F channel with L = 10 mkm 
and S* 5 mkni^ shown in Figure [1] with (a) n-doped and (b) p-doped silicon with 
high carrier concentration = 10^® cm""^. On each plot the three families of curves 
correspond to three positions of the "spin valve" emitter while the carriers at the 
collector are always highly polarized with a = 0.8. Each family has two curves labeled 
1 and 2 (bold) which corresponds to the low and high ratios of the junction and channel 
resistances ri/ro = 1/5 and ri/ro = 5/1. The solid curves Nl, N2 describe the case 
of non-polarized carriers at the emitter with /3 = 0.5 which was considered in the 
dashed curves PI, P2 are for the emitter polarization fi — 0.7 which is close to the 
collector polarization a = 0.8, and the dash-dotted curves A1,A2 are for the emitter 
polarization /3 = 0.3 which differs significantly from the one of the collector. For both 
type of carriers and for both low and high ratios of the contact/junction resistance 
the current saturation is observed being the indication of the spin blockade regime. 
The switch between different emitter polarizations leads to strong modulations of the 
current saturation amplitude, creating the possibility of the magnetization control on 
the current. 



least one ferromagnet polarization, say the emitter polarization /3. Below we shall see 
that these expectations are confirmed by the current-voltage dependencies in both n- 
doped and p-doped channels with both low- and high-Ohmic resistance. 



3. Current- voltage characteristics 



The current-voltage characteristics for n-doped and p-doped silicon channels with high 
and low carrier concentration are shown in Figures [3] and HI respectively. 
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The data for the room temperature drift mobihty of electrons and holes in the bulk 
silicon doped with various concentrations was taken from the standard reference plots 
[19] . allowing to estimate both channel resistance for a bulk silicon sample with L = 10 
mkm and S = 5 mkm^ and the reference critical current density from (fTTj) . In 
Figure [3] the results are shown for the silicon channel with high carrier concentration 
= 10^^ cm~^ where the reference critical current density j^^-* is within the range of 
160 . . . 180 mkA/mkm^, and the channel resistance r^ = 25 ... 30 Ohm. On each plot 
the three families of curves correspond to three positions of the "spin valve" emitter 
while the carriers at the collector are always highly polarized with a = 0.8. Each 
family has two curves labeled 1 and 2 (bold) which corresponds to the low and high 
ratios of the junction and channel resistances ti/tq = 1/5 and ri/rg = 5/1. The solid 
curves Nl, N2 describe the case of non-polarized carriers at the emitter with (3 = 0.5 
which was considered in [H], the dashed curves PI, P2 are for the emitter polarization 
P = 0.7 which is close to the collector polarization a = 0.8, and the dash-dotted curves 
Al,A2 are for the emitter polarization (3 = 0.3 which differs significantly from the one 
of the collector. For both type of carriers and for both low and high ratios of the 
contact /junction resistance ratio the current saturation is observed being the indication 
of the spin blockade regime. The switch between different emitter polarizations leads to 
strong modulations of the current saturation amplitude, creating the possibility of the 
magnetization control on the current. 

To compare the results for low-Ohmic samples with the high-Ohmic ones, we 
present in Figure H] the current- volt age characteristics for the silicon channel doped 
with low carrier concentration A^ = 10^*^ cm~^ where the reference critical current 
density is much lower and is within the range of 0.002 . . . 0.003 mkA/mkm^. The channel 
resistance is correspondingly much higher, and we consider tq = 200 . . . 500 kOhm for 
n-doped and p-doped samples, respectively. The labeling of all curves is the same as 
in Figure [3l and the maximum achievable currents are much lower due to significantly 
smaller carrier concentration which decreases the conductivity. One can draw the same 
conclusions about the spin blockade manifestation and polarization control of the current 
amplitudes mentioned above which obviously can be applied also for the case of low 
carrier concentrations. We can conclude that the main goal of the present model which 
is the achievement of a deep current modulation by the magnetization switch at the 
emitter can be reached in both low- and highly-doped silicon samples, with both n- 
type and p-type doping and with both low and high ratios of the junction/channel 
resistances. Hence, the proposed model of the tunable current-voltage characteristics in 
a ferromagnet-silicon-ferromagnet channel seems to be applicable to a rather wide range 
of ferromagnet-semiconductor structures. 

4. Conclusions 

We have studied the steady-state current-voltage characteristics of ferromagnet-silicon- 
ferromagnet channels with long bulk silicon sample having the length exceeding the 
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Figure 4. Same as in Figure[3]but for tlic (a) n-doped and (b) p-doped silicon channels 
with low carrier concentration N = 10^^ cm^'^. The labeling of curves is the same as in 
Figure m and the maximum achievable currents are lower due to significantly smaller 
carrier concentration which decreases the conductivity. The same conclusions about 
the spin blockade manifestation and polarization control of the current amplitudes 
mentioned in the caption for Figure [3] can be applied here for the case of low carrier 
concentrations. 



spin diffusion length. The current behaviour was investigated in the presence of spin 
blockade regime at the collector junction, and the dependence of the critical current on 
both collector and emitter polarizations has been obtained analytically. It was found 
that the current amplitude can be effectively tuned by varying the difference between the 
collector and emitter ferromagnet polarizations which allows to perform the magnetic 
manipulation on the electrical current in wide class of both n- and p-doped, low- and 
high-Ohmic semiconductor channels coupled to ferromagnetic leads. 
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